Stability of a granular layer on an inclined "fakir plane" 
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Abstract. - We present here experimental results on the effect of a forest of cylinder obstacles 
(nails) on the stability of a granular layer over a rough incline, in a so-called "fakir plane" con- 
figuration. The nail forest is found to increase the stability of the layer, the more for the densest 
array, and such an effect is recovered by a simple model taking into account the additional friction 
force exerted by the pillar forest onto the granular layer. 
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Introduction. — The flow and stability of granular 
layers on inclined planes have been widely studied in the 
last few years because it has a great interest from both 
ractical and fundamental point of view [T] . From a prac- 
tical point of view, there is for instance an important need 
in risk management of rocks and snow avalanches for a 
better prediction of the stability and flow rules of gran- 
ular media on rough inclines [2J[3]. From a fundamental 
point of view, the transition from the "solid" to the "liq- 
uid" state of granular media is still an open question, and 
the inclined plane is a widely useful tool for investigating 
the unknown rheology of such complex matter pQ . Numer- 
ous studies have been done concerning the stability of such 
a granular layer close to the flow transition [5T-H4J . A key 
result is that the stability of a layer of height h on a rough 
incline of slope is characterised by two limiting curves 
in the (h, 0) parameter space corresponding respectively 
to the spontaneous starting and stopping of the flow. Be- 
low the stopping curve, the granular layer is stable and no 
avalanche can develop. Above the starting curve, the layer 
is unstable and must flow. In between the two curves, the 
granular layer can be either in a "liquid" (flowing) or a 
"solid" (at rest) state, and avalanches can be triggered by 
small perturbations, that have to be larger when closer to 
the stopping curve [S]. The key point is that tan#, which 
corresponds to the ratio of the tangential to the normal 
forces, is related to the effective friction coefficient of the 
granular layer with the bottom plate. Thus exploring such 
stability and flow of a granular layer on an incline gives ac- 



cess to the effective granular friction forces [71TL3] . The two 
starting and stopping curves follow the same trend with 
a shift of only a few degrees in 0. For large h values, the 
two curves go asymptotically to the values corresponding 
to the two characteristic angles of stability of "infinite" 
granular piles: the maximum angle of stability and the 
angle of repose. The typical shape of the two curves are 
qualitatively the same whatever the experimental config- 
uration, such as the roughness of the plane (velvet clothe 
[5], sand paper [13], glued beads [M9]) or the flowing gran- 
ular material (spherical beads [3|6l|9l[T3] , sand [6l[T3], or 
even more anisotropic material such as copper [T3] ) . These 
starting and stopping curves have been recovered by dis- 
crete numerical simulations with collisional and frictional 
forces [HIE] and different theoretical approaches based 
on physical arguments [TTJtIB]. The theoretical prediction 
|17j for the link between the stopping height at a given 
angle 6 and the correlation length between flowing grains 
have been put in light both experimentally [ 1 OJ and with 
discrete numerical simulations [14j . In this Letter, we in- 
vestigate the influence of a forest of cylindrical obstacles, 
experimentally a forest of nails, on the stability of a gran- 
ular layer onto a so-called "fakir plane" . We present first 
experimental results and then a model that allows one 
to recover the essential of the experimental observations. 
The findings may be of importance in risk management 
for the possible stabilisation of snow avalanches or rocks 
for example by a forest of trees. 
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Fig. 1: (a) Sketch of the experimental setup and notations 
used, (b) Sketch of the fakir plane (top view). 



Experimental setup. — The experimental setup con- 
sists of a plane that can be inclined by an angle 9 from 
the horizontal up to about 45 degrees [fig. (Ha)]. The 
plane is covered by a velvet clothe which is chosen so 
that the grains have a larger friction with it than be- 
tween themselves. Nails of diameter D = 2 mm and length 
H = 35 mm are hammered in the plane (perpendicular and 
tips up) along a square lattice with a distance A between 
the nearest neighbours as shown in fig. QJ6), leading to 
a so-called "fakir plane" . Two different fakir planes have 
been built with A = 10"\/2 mm ~ 14 mm and A = 5%/2 mm 
~ 7 mm {A/D = 7 and 3.5). The particles used in the 
experiments are sieved glass beads (density p ~ 2500 kg 
m ) of diameter d g ranging from 0.13 mm to 1.2 mm with 
a relative dispersion in diameter of order 0.1. The cor- 
responding diameter values have been measured with a 
particle analyser (Morphologi G3, Malvern). The dimen- 
sions of the fakir planes used in the present study, with 
a length L = 300 mm and width I = 180 mm, are sig- 
nificantly smaller than the ones used in previous studies 
El E1H2JII3]- Nevertheless, we have checked this does 
not change significantly the results in the "non-fakir" case 
(A = oo) and believe that the results would not be af- 
fected significantly with larger dimensions in the "fakir" 
cases. It is worth noting that the granular layer deposit on 
the plane is never higher than the nails so that the height 
H of the nails can be considered as infinite and does not 
play any role here. Note also that the "fakir plane" just 
refers to the used plane hammered with nails but does not 
refer to any possible "fakir" state of materials laying at the 
top of nails rather than on the bottom of the plate as ex- 
hibited by liquids drops in an induced "hydrophobic" state 
by a forest of microscopic pillars [18] . The different exper- 
iments presented here correspond to the following ranges 
for the cylinder/grain size ratio D/d g and the number of 




1.5 < D/d g < 



15 



grains between two nails (A — D)/d g 
and 4 < (A - D)/d g < 100. 

The experimental procedure consists in pouring a 
known amount of grains onto the horizontal plane of sur- 
face S = L X I, and then slowly inclining it to trig succes- 



Fig. 2: Dimensionless layer thickness h/d g as a function of 
the angle 9 for the non-fakir plane and different grain sizes d g 
with open (solid) symbols for the starting (stopping) curves: 
(V,T) d g — 0.33 mm, (o, •) d g = 0.75 mm and (A, A) d g = 

1.2 mm. ( ) Best fits of starting data by eq. |l} with fitting 

parameter values (V) 6i m — 27.6°, 9 2m — 59.1° and A m ~ 
0.31mm, (o) 9 lm ~ 25.4°, 9 2m ~ 46.9° and A m ~ 0.71mm, 
(a) 9 lm ~ 19.1°, 9 2m =± 36.4° and A m ~ 1.1mm. (— ) Best 
fits of stopping data by eq. (JlJ with the fitting parameters: 
(T) 6i T ^ 23.1°, 9 2r ^ 58.8° and A r ~ 0.31mm, (•) 9 Xr ~ 
22.7°, 6 2r ^ 46.3° and A r ~ 0.79 mm, (a) 6 lr ~ 17.3°, 6 2r ^ 
35.5° and A r ~ 1.1mm. 



sive avalanches that affect the entire surface, following the 
procedure described in [S]. The weighting of the grains 
that have flowed out of the plane after each avalanche at 
an angle 9 allows one to infer the thickness h of the re- 
maining layer assuming a constant solid fraction <p of the 
grain layer (</> ~ 0.6). This procedure allows to reconstruct 
the two curves (h, 9) that characterise the stability of the 
granular layer on the inclined plane: the starting curve at 
which the layer starts flowing and the stopping curve at 
which the layer stops flowing. 

Experimental Results. 

Inclined plane without pillars. Let us first present our 
measurements for the "non- fakir" case (A = oo). Fig- 
ure [5] shows the starting and stopping curves in the (h, 
9) parameter space obtained for different grain diameters. 
These curves can be fitted by the following equation where 
the indexes m (for "motion" ) and r (for "repose" ) refer to 
the starting and stopping cases respectively: 



tan(9 mr (/i) = tan#i m . r + 



tan 9? 



tan( 



1 + hj \ m ,r 



(1) 



In this equation 9\ m , r corresponds to the characteristic 
starting or stopping angle of the granular pile (layer of in- 
finite height h — > oo), 92m, r corresponds to the characteris- 
tic starting or stopping angle for vanishing layer thickness 
h(h—>0), and A mjI . is the characteristic height of influence 
of the bottom rough condition for the starting or stopping 
case respectively. The indexes m, r will be dropped in 
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Fig. 3: Fitting parameters of eq. {TJ: (a) tan#i mjr (A, A) and 
tan 92m,r (°, •), and (b) \ m ,r as a function of the grain diameter 
d g for starting (open symbols) and stopping (solid symbols) 
cases. ( ) Guidelines for the eyes. ( — ) Linear fit X m ,r = d g . 

the following to avoid too heavy notation but one have to 
keep in mind that 6-y, 62 and A are not the same for the 
starting and stopping cases. Such an equation was already 
used by [TTfTS] for example, but an exponential dependence 
exp(— h/X) instead of the (l + h/X)" 1 dependence has also 
been used by [5ll6l[9l[T4] . Both functional dependences co- 
incide at low h/X values, giving a linear decrease of tan 6 
as a function of h, but differ in the asymptotic relaxation 
of 8 towards 8\ for large h values. The fits of our data by 
eq. H]) appearing in solid and dashed lines in fig. [2] for the 
stopping and starting curves respectively are quite close 
to the experimental data. The fitting values of the three 
parameters tan 8\ , tan 82 and A are shown in fig. [3] for all 
the used grain sizes d g from 0.13 to 1.2 mm. The tan#i 
value lies in the range 0.3 < tan0 x < 0.5 (17° < 6 X < 28°) 
with the starting value 0\ m above the stopping one 8\ r 
by about 2 degrees. No systematic variation is observed 
for tan 6*i with d g . By contrast, the tan# 2 value decreases 
significantly with increasing d g from about 2 down to 0.7 
(35° < 62 < 65°) with the starting value 02m above the 
stopping one 6*2 r by also about 2 degrees. This mono- 
tonic dependence of 82 found here for glass beads on vel- 
vet clothe is different from the non-monotonic dependence 
observed by [§] for glass beads on a plane with glued beads 
of a given size. The difference between tan #2 and tan#i 
here decreases which means that the effective friction co- 
efficient between the grains and the rough plane is less 
different than the effective friction coefficient between the 
grains when the grain size is larger. The parameter A m . r 
is found to increase linearly with the grain diameter d g 
with no significant difference between X m (starting case) 
and A r (stopping case) values. Note that the fitting value 
Xm.r/d g ~ 1 we find here for our case of glass beads can 
be significantly different for non-spherical particles (T3"] . 

Inclined plane with pillars. The influence of a forest of 
cylindrical obstacles (pillars) on the starting and stopping 
curves has been studied by performing similar experiments 
on the fakir planes with two nail spacings A and different 
grain sizes d g . Figure 0] shows the typical effect of the nail 
forest on the (a) starting and (6) stopping curves for grains 
of diameter d g = 0.75 mm. The (h, 8) curves obtained for 
the fakir planes present the same trend evolution as the 




Fig. 4: Dimensionless layer thickness h/d g for (a) starting and 
(b) stopping cases as a function of the angle 9 for grains of 
diameter d g = 0.75 mm and for (o, •) A = 00, (□, ■) A = 
14 mm and (0, ♦) A = 7 mm. ( — ) Best fits of starting data by 
eq. (TTJ with the fitting parameters: A m ~ d g , 62m — 46.9° and 
(o) 6? m ~ 25.4°, (□) 9 lm ~ 26.7°, (0) 9i m 31.1°. (— ) Best 
fits of stopping data by eq. ([TJ with fitting parameters values 
A r ~ d g , 6 2r ^ 46.3° and (•) 6£ ~ 22.7°, (■) lr ~ 24.8°, 
(♦) 9i r -29.8°. 

non-fakir plane. Nevertheless, for thick enough layers at 
small enough angles, a significant shift of the curves is 
observed: for a given layer thickness h, a higher 6 value 
is measured for the fakir case with a greater shift for the 
denser network (smaller spacing A) ; and for a given angle 
8, a thicker layer h is measured for a smaller A. This 
clearly shows the stabilising effect of the nail forest. This 
stabilising effect vanishes for vanishing layer thickness at 
large angles where the different curves collapse. In order 
to quantify the influence of the nail forest on the layer 
stability, we fit the different curves for the different planes 
with the same eq. ([1]) by keeping about constant the two 
parameters 82 and A that are related to the influence of 
the bottom wall and already extracted in the case without 
obstacles. The obtained fitting curves, shown in fig. H] 
by dashed (a) or solid (b) lines for the starting (a) or 
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Fig. 5: Evolution of tan#i relative to its value without pillars 
tan#f° (A = oo) as a function of the grain diameter d g for 
the two fakir planes with nail spacing A ~ 14 mm (□, ■) and 
7mm ((>,♦)■ Open and solid symbols are related respectively 
to the starting and stopping cases. ( — ) Linear fits with slopes 
0.07 mm -1 for A ~ 14 mm and 0.22 mm -1 for A ~ 7 mm. 

stopping (6) data, with the only one remaining free fitting 
parameter 0\ are quite good. The effect of the pillar forest 
on the layer stability thus reduces to the modification of 
the Q\ value. The dependence of tan0i for the starting 
and stopping cases relative to its value tan0j > ° without 
any pillars (at A = oo) is plotted in fig. [5] for the different 
grain diameters d g and the two pillar networks A. One can 
see that the difference tan0i(A) — tanf?^ is proportional 
to d g with a slope that increases for decreasing A, with no 
significant differences between the starting (open symbols) 
and stopping (solid symbols) cases. 

Modelling. — Let us now model the stabilisation 
played by the forest of pillars. It is known that a granu- 
lar flow around one cylinder of diameter D induces a drag 
force F c on this cylinder [HOUD]. For a quasi-stationary 
regime, the force exerted by the cylinder forest on the 
granular flow should be opposite to the force exerted by 
the granular flow on the cylinder forest, and we shall con- 
sider it to be equal to the one-cylinder force times the 
number of cylinders by neglecting any strong flow inter- 
action between the cylinders. The stress field of the grain 
contact network around a cylinder in relative motion has 
not been yet investigated, but the granular flow pertur- 
bation has been shown to be localised in a radial exten- 
sion zone of about one cylinder diameter whatever the 
grain size [50]. The one-cylinder force F c is known to 
be F c = aptfigDh 2 independent of velocity, where p<j) is 
the effective density of the granular medium (0 refers to 
the grain packing fraction), h s is the vertical penetration 
depth of the cylinder from the granular free surface, and 
a is a dimensionless coefficient depending on geometrical 
shape, and on friction and restitution coefficients of the 
grain/cylinder couple [TO]. Note that some weak depen- 
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Fig. 6: Same data as in hg. [5] but as a function of the ratio 
Ddg/A 2 . ( — ) Linear fit with slope 5.5. 

dence of F c with the grain diameter d g (or the size ratio 
dg/D) has been observed by some authors (THIHO], which 
will here be neglected in a first approximation. As for 
a forest of cylinders in a square array of spacing A the 
number of cylinders per unit area is 1/A 2 , the force per 
unit surface exerted by the cylinder forest onto the gran- 
ular layer is F c f = F c /A 2 = apcf>gDh 2 /A 2 . If one now 
considers the stability of a layer of thickness h s close to 
the surface of a granular pile and far away from any bot- 
tom, the stability of such a layer of length L and width I 
when it just starts or stops at the characteristic starting 
or stopping angle 6\ can thus be written as 

= p(f>gh s lLsm0i — fi g p(j)gh s lL cos B\ 

— —^apifigDh 2 cos 6\, 
A z 

where p g = tan^J corresponds to the internal friction of 
the granular packing and thus to the characteristic angle of 
the pile for the starting or stopping avalanche flow referred 
usually as the "maximum angle of stability" (0^) and 
"angle of repose" of the pile. Equation ^ leads 

thus simply to 

tan 6»i (A) = tan 0f° +a-r^. (3) 
A z 

The additional friction term [second term in the right hand 
side of eq. fl3J)] coming from the cylinder forest increases 
thus the characteristic starting and stopping pile angles 

01m, r- 

If one now considers the configuration of a granular layer 
of thickness ft, on a fakir plane, its stability can thus be 
written as 

» n / a \ tan 09 ^ tan 0i (A) . , 

tan 9(h, A) = tan 0i (A) + , v ' . (4 

1 + h/X 

This equation is the same as eq. ([T]) but with the effect 
of cylinder forest encoded here in the tan0i term which 
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is a function of the cylinder spacing A (and diameter D) 
given by eq. ©. The experimental results plotted in fig. [5] 
shows that tan#i — tan0J° scales with d g and thus the 
additional friction term in eq. ([3]) scales with d g . This 
means that the characteristic flowing thickness h s consid- 
ered in model equation (3) scales with d g and docs not 
correspond to the entire layer height h. By considering 
thus that h s oc d g , model equation (j3J predicts that the 
difference tan#i — tan(9J° is proportional to the size ra- 
tio Dd g /A 2 . The plot of our experimental data corre- 
sponding to tan 0\ — tan#J° as a function of Dd g /A 2 in 
fig- El appears about linear with the slope 5.5. This slope 
value gives thus ah s ~ 5.5<i g which gives in turn h s ~ d g 
when considering the value a ~ 5 ± 1 measured in |19) . 
The fact that the found characteristic flowing thickness 
h s scales with the grain diameter d g is not very surprising 
as granular matter is known to exhibit complex fluid be- 
haviour with flow localisation [TJ. In particular flows are 
very localised in a thin surface layer as already visualised 
by [21] for unsteady "spontaneous" avalanches or by [22] 
and references herein for forced flows: in the first case, 
the velocity profile exhibits an exponential decrease in the 
depth with a characteristic flowing layer that corresponds 
to about (2± l)d g [5T| ; in the second case, the linear part 
of the velocity profile related to the high granular flow is 
followed by an exponential decreasing part with a charac- 
teristic thickness (1.1 ± 0.1)<i g [55]. Besides, the typical 
flowing thickness deduced from the modelling of the pile 
stabilisation by close lateral walls have also been shown to 
be of about (2 ± l)d g for the range of grain diameter used 
in the present study where the grain size is not too small 
to be hardly influenced by van der Waals interaction forces 
[23,23]. It is worth noting that h s is found here similar 
to A, as A ~ d g and h s ~ d g . This means that the effec- 
tive friction and thus the stability of a granular layer on 
a rough plane may be related to the avalanche flow locali- 
sation. The granular flow may reach the bottom plane for 
low enough h and not for high enough h. This may explain 
the shape of the h(6) curves of figs. 2 and 4 where tan0 
varies strongly at low h/d g and weakly at large h/d g . Note 
that this would imply that the angle of the surface granu- 
lar layer 8 S may be different from the plane angle for large 
enough granular layer h/d g as the granular flow does not 
"feel" the bottom anymore: the granular layer would be 
no more homogeneous along the plane but would exhibit 
a streamwise negative gradient dh/dx = — tan(0 s — 8). 
In the present study, we restrict to h/d g < 20, where no 
significant thickness gradient is observed along the plane. 
For very large h/d g values (h/d g > 20), we indeed observe 
such a thickness gradient. In such a case, we believe that 
eq. (4) still holds but with the angle of the surface layer 
9 S instead of the plane angle 8. Such a regime have al- 
ready be seen between two close lateral walls for which 
the pile angle is stabilised by the additional friction forces 
at the walls [53] and may lead to the so-called "super sta- 
ble heap" for a large enough forced flow [53] . The effect of 
the pillar forest on the pile stabilisation appears to be very 



similar to the effect of close lateral walls, both cases lead- 
ing to an effective additional frictional term. Note that in 
the explored range of parameters, we always observed an 
angle 0\ smaller than 82 but the contrary may be imag- 
ined for a high enough value of the additional friction term 
due to the pillar forest in eq. (3). Note also that when h 
is smaller than h s , the additional stabilising term of the 
cylinder forest in eq. (3) is thus overestimated, but this 
term is negligible as the stabilising effect of the bottom 
is much larger and governs the layer stability in such a 
case. Moreover, since the measured drag force on a cylin- 
der has been reported by [ini[50] to increase weakly with 
decreasing grain size d g (or with the size ratio d g /D), the 
variation of tan#i(A) — tan#i with d g and also Dd g /A 2 
would be perhaps weakly sublinear. 

The critical granular thickness necessary to observe a 
significant increase of its stability by a given pillar network 
can be estimated by the cross-over between the regime 
dominated by the bottom friction given by the second term 
of the right hand side of eq. (U) and the regime dominated 
by pillar forest friction given by the second term of the 
right hand side of eq. ([3]). As the variations of tan#i 
are small compared to the difference tan 82 — tan 8\ , this 
critical height value h c is given by 

tan 82 — tan 8\ Dh s 



As this happens typically when h/X > 1, and as A ~ d g 
and h s ~ d g , the previous equation leads approximately 
to 

h c tan 9 2 — tan 9% A 2 

1 7m~' W 

d g a Dd g 

This predicts that the critical dimensionlcss layer thick- 
ness h c /d g scales linearly with the dimensionlcss size ra- 
tio A 2 /(Ddg) for a constant ratio (tan 82 — tan#i)/a, and 
thus that the critical layer thickness h c scales linearly with 
A 2 /D. Note that (tan 82 — tan#i)/a should not be con- 
stant but vary only weakly with d g as both tan 82 — tan 8\ 
and a decrease with increasing d g . 

In this paper, we show a clear stabilising effect of an ar- 
ray of pillars for a granular layer on a rough incline. The 
experimental results are well fitted by a model equation 
taking into account the additional friction arising from the 
pillar forest when compared to the usual bottom friction 
on the rough incline. The key parameter appears to be 
the dimensionless size ratio A 2 /(Dd g ). Additional exper- 
iments would be necessary to verify the dependence on 
the pillar diameter D and the model prediction that the 
real flowing thickness h s should be smaller than the layer 
thickness h for high enough h and should scale with the 
grain diameter d g . In forthcoming experiments, we in- 
tend to look at continuous feeding flows where stationary 
regimes may be reached. Finally, the present study may 
be useful for risk management as it predicts that a forest 
of trees characterised by a typical spacing A and diameter 
D may stabilise some layer of discrete materials such as 
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rocks or even snow on a slope. It is also related to the 
possible stabilisation of soil slopes by vegetation and its 
root system. 
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